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Abstract 

We study the possibility of localization of the front present in a one-dimensional 
ballistically-contr oiled annihilation model in which the two annihilating species are 
initially spatially separated. We construct two different classes of initial conditions, 
for which the front remains localized. 
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1 Introduction 



During the last decade, a large body of work has been devoted to the study of 
the kinetics of diffusion-annihilation processes. It is now well established that, 
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below some upper critical dimension, the fluctuations play a central role and 
that accordingly mean-field like approximations are inappropriate. 

Moreover, for the two species case A + B — > 0, when the two reactants are 
initially spatially separated, a reaction-diffusion front, getting larger with time, 
is formed [1,2]. In the long time regime, the time dependent properties of this 
front (position, width) are characterized by power laws which are non-mean- 
field at or below two dimensions. 

Instead of investigating a time dependent problem, it was shown by Cornell 
and Droz [3] that it may be advantageous to study the front formed in the 
steady state reached by imposing antiparallel currents J a = \ J\ and Jb = —\J\ 
of A- and B-particles at x — — oo and x = +oo respectively. It turns out that 
first, exact prediction can be made in the one dimensional case and second 
this stationary problem is closely related to the time- dependent one. 

More recently, a different but related problem has been considered, namely 
the case of ballistically-controlled annihilation processes. Most of the results 
obtained are for the one dimensional case. Initially, the particles are randomly 
distributed in space and their velocities are distributed according to a given 
distribution. The particles move freely and when two of them collide they 
instantaneously annihilate each other and disappear out of the system [4-8,12]. 
In 1985, Elskens and Frisch [4] considered the case of one species of particles 
moving with velocities +c and — c. The case of an arbitrary discrete velocity 
distribution has been studied by Droz, Rey, Frachebourg and Piasecki [7,8] and 
exactly solved for a symmetric three velocities distribution, using an exact 
closure of the hierarchy obtained previously by Piasecki [6]. Such processes 
can model several physical situations as a recombination reaction in the gas 
phase or the fluorescence of laser excited gas atoms with quenching on contact 
(the one-dimensional aspect can be obtained by working in a suitable porous 
media [9]) or, the annihilation of kink-antikink pairs in solid state physics [10]. 

We have recently investigated the case of the front formation for ballistic 
annihilation with two species initially separated in space in one-dimension [12]. 
This process can model, for example, the situation in which chemical species 
incorporated in a gel move ballistically under the action of a drift [11]. As 
the two species cannot penetrate one into the other (as they annihilate on 
contact), a well defined reaction front is formed. For a symmetrical Poissonian- 
type initial spatial distribution, we have proved [12] that the front behaves like 
a random walker for long time. In other words, at long time, the front can be 
anywhere on the line: it is delocalized. 

In view of what has been done for the diffusive case, it is natural to try to find 
a system with ballistic-annihilation possessing a localized reaction front. Let 
us first remark that in the ballistic case, to impose a flux of particles at the 



2 



boundaries of a finite system is equivalent to fix an initial spatial distribution 
of a large system. For example, the case with homogeneous spatial Poisson 
distribution can be mapped onto a system with boundary fluxes of particles 
at constant rate, the time between two particles input being exponentially 
distributed. 

Thus we shall consider the problem in terms of initial distribution of particles 
and ask the question: is it possible to localize the front by choosing a suitable 
initial distribution of the particles ? The answer is obviously yes. Indeed, as a 
simple example, consider the case in which initially each particle is located on 
a regular lattice, the A particles being to the left of the origin and the -B's to 
the right, in a symmetric way. The position at which collision between one A 
and one B particle takes place defines the position of the annihilation front. 
Thus the front will always sit at the origin. However, this case is pathological in 
the sense that the initial state is totaly ordered. For more general situations, 
with spatial randomness in the initial state, we have first to define what is 
meant by "localisation of the front". Then we shall show that indeed, it is 
possible to find suitable initial conditions localizing the front. 

The paper is organized as follows. In section 2, we define precisely the class 
of models studied and introduce several criterions of localisation. In section 
3, we consider a class of correlated initial conditions for which we are able 
to compute the probability density /i(X] t) to find the front at position X at 
time t and which leads to the localization of the front. In section four, we 
consider a different class of initial conditions with no correlations between the 
initial positions of particles. For this more physical class of initial conditions, 
we are not able to compute explicitly p,(X] t), but we can show that the mean 
square position of the front (X 2 ) converges for t — * oo. This means that the 
probability density p(X; i) does not spread out. Remarks and conclusions are 
drawn in section 5. 



2 The model and the criterion of localization. 

We consider a one-dimensional system formed of two species of particles. Ini- 
tially, particles A are spatially randomly distributed in the region (— oo, 0) and 
the B ones are spatially randomly distributed in the region (0,oo). For the 
sake of simplicity we shall consider only symmetric distributions in A and B. 
In particular, pa(— x) = Pb{x) = p(x),x > 0, where pa and pn are the number 
densities of particles A and B respectively. For nonsymmetric distributions, 
the front moves acquiring some mean velocity. 

The velocity of each particle is an independent random variable taking the 
value ±c with equal probability. Particles of the same kind suffer elastic col- 
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lisions. When two particles A and B meet, they annihilate. Thus, practically 
the A particles with velocity — c and the B particles with velocity +c move 
freely. Accordingly, the relevant part of the dynamics concerns only the A 
particles with velocity +c and the B ones with velocity — c. 

Let y2, ■ ■ ■ , Vk, ■ ■ •) be the initial positions of the A particles with velocity 
+c and (xi,X2, ■ ■ ■ ,x k , ■ ■ .) the initial positions of the -B's with velocity — c. 
The particles are labeled such that yk < • • • < y\ < < x\ < ■ ■ ■ < x k - The 
relative velocity between the A and B particles being 2c, the pair of particles 
initially at {xk,yk) will collide at time tk = (xk — y k )/2c. The position at which 
this collision will take place defines the position of the annihilation front. Thus, 
the position of the front at time t is (xi — yi)/2c, if let < x± + yi (i.e. if no 
particle has yet collided), (X2 — ?/2)/2c, if x\ + y± < let < x 2 + t/2 (i-e. only 
the first particles on the right and on the left have collided), and so on. The 
dynamics in itself is purely deterministic, the only stochastic aspect comes 
from the initial conditions. 

The properties of the front are completely defined by the probability density 
/i(X;t) to find the front at the point X, at time t. In a previous paper [12], 
we have proved that: 

n(X;t) = {5(X-±[ Xl + yi ])) (1) 

+ (jt [ 6 ( X - 2^+1 + Vk+i\) ~ 6 ( x ~ \W + Vk])\ 0(2ct - [x k - y k })^ , 

where S and 9 are the usual Dirac and Heaviside functions and the brackets de- 
note the average over the initial positions, according to the initial distribution. 
It follows that, the second moment of yu(X; t) is 




xi + yi 



^ x k +i +l/fc+i x 2 



x k + Vk 



9 (2ct - [x k - y k \) 



(2) 



Several criterions can be introduced to define the localization of the front. The 
first is simply to ask that fi(X; t) approches for t — > 00 a limiting distribution 
with finite moments. A second criterion consists in asking that /i(X] t) has 
a compact support with respect to X, for all time. These two criterions are 
useful provided that one is able to compute /j,(X;t). This calculation will be 
possible for our first class of initial conditions (section 3). 

However to deal with the second class of initial conditions, we are led to 
propose a weaker definition of the localization. The front will be said to be 
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localized if fi(X; t) has a finite second moment (or variance in case of a non- 
symmetric initial distribution) for all time t < oo. 

We know what is happening in two limit cases. For a Poisson distribution the 
front is not localized while for particles regularly distributed on a lattice the 
front is strongly localized. It is thus natural to consider an initial distribution 
which extarpolates smoothly between these two limits. Such distribution is 
provided by the Erlang-k-process [13]. For k — 1, one recovers the Poisson 
process while for k — > oo one finds the case of the regular distribution. How- 
ever, it can be shown [14] that for all finite values of k the front is delocalized. 
Only the limit case k — > oo is localized. 



3 Correlated initial distribution. 



We aim at finding a class of initial conditions for which the probability den- 
sity /i(X] t) has a compact support, with respect to X. This property can be 
obtained with the following initial condition: each B particle is uniformly dis- 
tributed in an interval of length 5 < £, centered at (k — where £ is the 
distance (center-to-center) between two consecutive intervals and k labels the 
particles (k = 1,2, . . .). We use the same distribution for particles A, but on 
the left of the origin. First, remark that the position of any particle is com- 
pletely independent of the other; in addition we choose identical distribution 
for each particles, but at different regularly spaced positions. We eventually 
conclude that /i(X;t) must be a time periodic function of period £/c. Hence 
it is enough to compute /i(X;t) for t G [0,£/c). In this case, eq. (1) simply 
reduces to 



fi(X;t) = u (X) + u(X;t) 

= (S(X-l[x 1 + y l })) (3) 
+ ([S(X -\[x 2 + y 2 }) - 5(X -\[x l + yi })] 6(2ct - [x, - Vl ])) . 

The first bracket reads 



i+S 
2 



e+s 

2 



e-6 

2 



4(f-mW-m). 



(4) 



The computation of the second bracket is less straightforward, mainly due to 
the presence of the Heaviside function: 
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M+S 



e+6 



u(X;t) = - J dy 2 J dx 2 J d Vl J dx 1 9(2ct-[x 1 -y 1 }) 

Sl+6 Zt-5 t+5 1-5 

2 2 2 2 

x [5(X -\[x 2 + y 2 \) - 5(X - \[ Xl + yi \) 



(5) 



However, simplifications arise by doing a Laplace-transform. Carrying the four 
integrations over y ly y 2 , X \ and x 2 and inverting, we finally obtain, v[X\ t) as 
a product of the Heaviside function 6(^5 — \X\) and a sum of five terms. 
We shall not reproduce these terms here as they are a little bit cumbersome. 
Nonetheless, we plotted both and v{X\ t) for t e [0,£/c) (see figure 1 

and 2). The plot of u(X;t) shows that it is not a simple function: already at 
this stage, t) has a non-trivial shape. 




Fig. 1. Plot of n(X; t), for £ = ±, 5 = \, X G [0, \] and ct € [0, 1]. The probability 
density /J-(X; t) is symmetric with respect to X and periodic (with period 1) with 
respect to ct. 



It is clear that if we take an arbitrary distribution instead of a uniform one 
for the probability density inside each interval, we should always find that the 
support of fi(X;t) is [— 1<5, |<5], although the detailed shape of fi(X;t) may 
vary considerably. Moreover, one can easily generalize all these statements to 
non-symmetric distributions, by taking into account that the front will move. 
The support of fi will not be the same for all time, as it will follow the mean 
position of the front. 
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Fig. 2. Plot of u(X;t), for t = \, S = \, X G [-±, \] and ct G [|, §]. 
4 Noncorrelated initial conditions. 



In this section we show that the front can be localized by a generalized Poisson- 
type distribution. For a homogeneous symmetric Poissonian particle distribu- 
tion, we found [12]: 

(X 2 ) ~ ^ + Const (t -> oo) , (6) 
2p 

where p is the initial density. Now suppose that instead of having a constant 
p, the density increases with position according to a power law: 

p = p(x) = x a . (7) 



Then, the front (whose average position is the origin) will see an increasing 
density in time. We can expect to obtain 

(X 2 ) ~ + Const (t -> oo) . (8) 



For a > 1, (X 2 ) converges to a constant, leading to a localization of the front, 
whereas for a < 1, it diverges and the front is delocalized. The mechanism 
leading to the this result is clear: when we increase the density, the mean 
distance between two successive particles decreases and thus two consecutive 
collisions should happen relatively often. As a consequence, as the time grows, 
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the probability to find the front far from its mean value becomes smaller. We 
shall now make more precise this qualitative picture. 



To begin we choose, as indicated, the initial particle distribution to be an 
inhomogeneous Poisson law. In other words, if there is a B particle at position 
Xi, the probability to find its right nearest neighbor at position x 2 > X\ is 



tt(x 2 ,x 1 ) = p(x 2 )exp 



X 2 

J p(z) dz 



(9) 



where p(x) should be understood as the initial particle density at point x. For 
p(x) = Const, we find again the homogeneous Poisson law already studied. 
The probability to find k — 1 particles between and x k and a particle at x k 



is 



ir k (x k )=p(x k ) J dx fe _i exp - J dz p(z) 

Sk_i 

Zfc-i Zfc-i 

xp(x fe _i) J dx k ^ 2 exp - J dzp(z 



Xk-2 



X2 %2 

x p(x 2 ) J dxi exp — J dz p(z) 



XI 



p(xi) exp 



XI 



dz p(z) 



L 



■p(x k ) exp 



dz p(z) 



Uo k dzp(z)] 
(k-l)\ 



k-i 



(10) 



n k is normalized to 1 provided the density p(x) does not vanish at infinity, 
which is satisfied for the law (7), with a > 0. The distribution of particles A 
is chosen in a symmetric way Pa(u) = p(~y), V < 0. 

Using eq. (2), the second moment becomes 



xi \ o 

( x2 )t = f dx 1 p(x 1 )exp i- J p\ J dyip(yi)exp 

/ -oo 



xi + yi 



yi 

Vk 



J dx k p(x k ) J dy k p(y k ) J dx k+1 p{x k+1 ) j dy k+1 p(y k+1 ) 

— oo 

Cfc+l 

/-/ 



+ 

k=1 -oo 

x exp | - y p - / p 

Vk+i 



[(*-!)!] 
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X 



fx k+ i + y k+1 \ 2 fx k + y k 



V 2 J 
where we used the notation 



9(2ct- [x k -y k ]) , 



J P = J dzp(z) . 







(11) 



Now integrating twice by parts the expression 



/ 



dx 







k+l 



dx 



k+l 



X 



Vk 

J dy k+1 







Xk+l 

<>X P \ j p 





dy k 



+i 



cxp I - j p 

Vk+l 



S x k+1 + y k+ i 
V 2 



and inserting the result into eq. (11), we obtain 



x y 

(X' 2 ) t = I dxp(x) I dyp(y)exp \ - J p- J p 

o o 





oo oo 



+ t[dxp(x)[dyp(y) { -Il^ 



x 





oo 



oo oo 

\j dx' j dy' 



We go on by choosing: 



J dx' exp I - j p - j p | (x' - y) 



fx' y 



cxp 



p- J p 

\ 



p{x) = x a , 



x + y 



9(2ct- [x + y]) . 



(12) 



(13) 



with a > 0. As long as the time is finite, {X 2 ) t is finite for any a, because the 
a = case appears to be an upper bound. However it can grow arbitrarily 
when t goes to infinity. Can we find a > such that (X 2 ) t is bounded for any 
time? It is very difficult to calculate exactly (X 2 ) t . Even when t is large, we 
have not been able to find its asymptotic expression. Nevertheless we can put 
t — > oo in the expression of (X 2 ) t and check its convergence. Clearly, as the 
first term of eq. (12) is time independent, we only have to care of the infinite 
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sum. In addition, when t — > oo, assuming that we can interchange the limit 
with the sum and the integrals, the restriction on the domain of integration 
introduced by the Heaviside function is removed. We are thus left to calculate 



x a y a / x a+l y a+l \ * 



a k = J dx J dy 



X 





oo 



oo oo 



X X 



(k\) 2 \a + la + l / 



' x' a+1 


y a+l 


v a + 1 


a + 1 


x' a+1 


y' a+1 N 


a + 1 


a + l, 



and to check the convergence of YJk=o a k- 
By changing the variables 



x 



a+l 



y 



a+l 



a + 1 ' a + 1 ' 



x 



/a+l 



a + l 



./a+l 



a + l 



(14) 



we get 



oo oo 



afc = / da / db 

o o 



(abY 



da' [(a + l)a']"^T e 



(k\) 2 

x {[(a + l)a']~^rr - [(a + l)6]~s+r} 

oo oo 

--/da'/di/ 



-o'-6 



^a + lya'b 



a+l g-a'-b' 



(15) 



We finally obtain (by changing the order of integration over a and a' and over 
b and 6'): 



a fc = (a + 1) 



r(fc + 2) 



A; + 1 + 



1 



2(a + 1) 



r 2 (aTT + fc + 1 

T 2 (A; + 2) 



(16) 



To study the convergence of Y^kLo a k, we simply have to study the asymptotic 
behavior of a k . We find (see [15]) 



1 , _2d 2_ 

a k = -k 



l--}— + —^—; + 0(k~ 1 ) 

a + l (a + l) 2 v 
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when k — > oo. It is thus evident that the sum converges if and only if 
Ik^M dk k~ 2+ ~ does (where M is any integer sufficiently large). This leads 
us to the previously announced result: (X 2 ) t converges if a > 1, diverges if 
a < 1. For a = 1, one has a logarithmic divergence. 



5 Conclusion 

We have exhibited two different classes of initial conditions leading to a lo- 
calization of the reaction front for this ballistic annihilation model. We have 
shown that the main physical reason for derealization is a too slow decay of 
the probability to find two consecutive particles at an arbitrarily large dis- 
tance. Even an exponential law characterizing the Poisson distribution is not 
enough. 

Several questions remain to be adressed. One natural extension of this work 
is to know whether or not a faster decay of the probability distribution of the 
nearest neighbours (for example a Gaussian probability instead of an exponen- 
tial one) would lead to localization. Indeed, we know that if the interparticle 
distribution is 

_I 2_I 2 
Tr(x 2 ,X 1 ) = X 2 e 2 X 2 2 x i 

(see eq. (9)), then (X 2 ) t diverges logarithmically. 
Perhaps, the distribution 

7r(x 2 ,a;i)=e-5(«-«) a > 

could lead to a front localization. Another possible extension is related to 
the case of a more general distribution of the velocities. These more difficult 
problems are under investigation. 
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